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Abstract. Let [M, g) be an odd-dimensional incomplete compact Riemannian singular 
space with a simple edge singularity. We study the analytic torsion on M, and in par- 
ticular consider how it depends on the metric g. If g is an admissible edge metric, we prove 
that the torsion zeta function is holomorphic near s = 0, hence the torsion is well-defined, 
but possibly depends on g. In general dimensions, we prove that the torsion depends only 
on the asymptotic structure of g near the singular stratum of M; when the dimension of the 
edge is odd, we prove that the analytic torsion is independent of the choice of admissible 
edge metric. The main tool is the construction, via the methodology of geometric microlocal 
analysis, of the heat kernel for the Friedrichs extension of the Hodge Laplacian in all degrees. 
In this way we obtain detailed asymptotics of this heat kernel and its trace. 



1. Introduction 

One of the key achievements in modern spectral geometry is the proof by Cheeger and 
Muller of the Ray-Singer conjecture, which equates the analytic and Reidemeister torsions of 
a compact smooth odd-dimensional manifold (equipped with a flat Hermitian vector bundle). 
Since one of these quantities is analytic and the other combinatorial, their equality has many 
important applications in fields ranging from topology and number theory to mathematical 
physics. The original definition of analytic torsion, and its conjectured relationship with Rei- 
demeister torsion, appeared in the famous 1971 paper of Ray and Singer |26j . The original 
proofs by Cheeger [11] and Muller [25] are still of great interest, but there are now several 
other proofs of this result as well, each with its own significance and leading to further gen- 
eralizations. Amongst these we mention in particular the ones based on Witten deformation 
[32] and 'analytic surgery' [18] . 

Singular spaces arise naturally in many parts of mathematics, and the development of 
analytic techniques to study partial differential equations on them is a central challenge in 
modern geometric analysis. Important examples of singular spaces include algebraic varieties 
and various moduli spaces, and they also appear naturally as compactifications of smooth 
spaces or as limits of families of smooth spaces under controlled degeneration. 

A natural and still-outstanding open problem is to determine whether the Cheeger-Miiller 
theorem has any analogue for compact stratified pseudomanifolds (which are perhaps the best 
behaved type of singular space). We refer to §2 of [2] for a detailed explanation of the differ- 
ential topological structure of these spaces. There are many difficulties even to formulate a 
precise conjecture. The natural homology theories in this setting are the intersection homol- 
ogy spaces of Goresky and MacPherson; however, there are many of these, none necessarily 
preferred over the others, which leads to an ambiguity in what one should mean by intersec- 
tion Reidemeister torsion. We refer to recent work by Dai and Huang [12J for a study of this 
issue for spaces with isolated conic singularities, which is already not so straightforward. 

It is not so easy to define analytic torsion in this setting either. The development of 
techniques to study elliptic differential operators, in particular geometrically natural ones 
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such as Dirac- or Laplace-type operators, on singular spaces goes back many years, starting 
from the work of Kondratiev and his school in the early 1960's. Cheeger [9], [10] was the first 
to understand the tractability of studying spectral geometry for smoothly stratified spaces 
endowed with 'incomplete iterated edge metrics'. This was the beginning of many further 
developments by various authors. The analysis of partial differential operators for spaces 
with conic (and the closely related asymptotically cylindrical) metrics was also studied by 
Lesch [20], Melrose [22], Briining and Seeley [6], Gil and Mendoza [15] and Mooers [21], to 
name just a few. Extensions to spaces with simple edge singularities were developed by the 
first author [21j . and Schulze and his collaborators [27j . see also [8]. Further extensions to 
spaces with iterated edge singularities are still in a less refined state of development, though 
see [2]. 

The present paper focuses on the definition of analytic torsion on the next simplest class of 
spaces beyond the ones with isolated conic singularities. Namely, we consider spaces (M,g) 
with a simple edge singularity and incomplete edge metric. This means the following. We 
assume that M is a compact stratified space, with a single top-dimensional stratum and only 
one other lower dimensional stratum B, which is therefore a smooth closed manifold. It is 
convenient, and not too misleading, to refer to the main stratum as M, and to write M 
when we wish to emphasize that we are talking about the entire space. We set b = dimB 
and m = dimM. The stratification hypothesis provides a neighbourhood IA of B which is 
the total space of a smooth bundle over B with fibre C(F), an open truncated cone over a 
compact smooth manifold F of dimension /, i.e. C(F) = [0, 1) x Fj ~ with (0,2:) ~ (0,2:') 
for any z,z' £ F. The metric g on M is arbitrary away from U, but in this neighbourhood 
takes the special form 

(1.1) g := dx 2 +x 2 k(x) +<j)*t)(x), 

where x is a smooth function on U \ B which restricts to a radial function on each conic 
fiber, <ft : {x = const.} — > B is the fibration of each level set, and n(x) is a family of smooth 
metrics on F depending smoothly on the parameter x for < x < 1, and t)(x) is a family of 
smooth metrics on B which likewise depends smoothly on x £ [0, 1). We give a more invariant 
definition of this class of metrics below. For simplicity we call (M, g) a simple edge space. 

Our goal is to study the analytic torsion of simple edge spaces. Ideally, one hopes to prove 
that the analytic torsion is well-defined and independent of the metric g, and that it defines 
an invariant which can be computed in terms of combinatorial (and perhaps other) data. We 
accomplish the first part of this here: namely we prove the existence of the analytic torsion 
for simple edge spaces and show its invariance properties under various conditions. 

To state our main result, let us recall some definitions and standard terminology. Let A^ ^ 
be the Hodge Laplace operator on a compact smooth manifold M with respect to the metric 
<7, acting on A:-forms, and det' its zeta-regularized determinant, with the zero modes removed. 
Define the (scalar) analytic torsion T(M, g) by 

m 

log T(M,g) = ^(-l) fc Hogdet'A fe)S . 

fc=0 

The determinant line on M is defined in terms of the de Rham cohomology by 

m (— l) fc 

(1.2) det^*(M):=(g)(/\ OP ^(M)) 

fc=0 

Identifying H^(M) = ker Afc i9 , then the determinant line on M inherits a natural L 2 Her- 
mitian structure || • ||^2 induced from the inclusion ker A^ 5 C L 2 Q, k (M, dV g ). There is another 



Rafe Mazzeo and Boris Vertman 



3 



Hermitian structure on this same line bundle, described by the norm 

II -11(2,3) -=T(M,g)\\ -\\ L2 ; 

this is called either the analytic torsion, the Ray-Singer or the Quillen metric. 

Now return to the setting where M has a simple edge. We shall introduce in §2 various 
classes of metrics on M. The most restricted is the class of rigid (incomplete edge) metrics, 
which are exactly warped product conic on each conic fibre of the tubular neighbourhood 
14, and amongst these the 'admissible' metrics which are rigid and are also a Riemannian 
submersion in IA. The more general types of metrics we study are asymptotically of this form, 
i.e. g = go + h where go is admissible and h is smooth (or poly homogeneous) up to the edge 
and decays relative to go- Finally, amongst these asymptotically admissible metrics we single 
out two special subclasses: 

i) The metric g is strongly asymptotically admissible if \h\ go = 0(x b+1 ); 

ii) The asymptotically admissible metric g is called even if the expansion of h (as de- 
scribed below) contains only even powers of x. 

The precise definitions are given at the end of §2.1. 
We may now state our main 

Theorem 1.1. Let (M,g) be a compact simple edge space. Then the analytic torsion, and 
hence the Ray-Singer metric \\ ■ \\ g)' ^ n ^ erms °f ^ e Friedrichs extension of the Hodge 
Laplacian, is well-defined for any asymptotically admissible metric provided m = dimM is 
odd, and moreover, \\ • ||(^ g ) is invariant under all deformations amongst strongly asymptot- 
ically admissible metrics g which fix the admissible rigid metric go- Ifb= dimB is also odd, 
then || • ||(|f s ) is invariant under all deformations amongst even or strongly asymptotically 
admissible metrics g, including those which vary the admissible rigid metric go- 

The earliest study of analytic torsion and its possible relationship to combinatorial invari- 
ants in the setting of singular spaces was by Dar [13] . More recently, the second author studied 
the analytic torsion of truncated cones [31], with similar independent work by Spreafico |30j. 
Dai and Huang [12] have initiated a study of Reidemeister torsion for conic spaces, but find 
that there is no obvious unique generalization of Reidemeister torsion. The theorem above 
suggests that this may be no accident, since the analytic torsion is invariant under a reason- 
ably broad class of metrics only when dimi? is odd, which excludes the case when B is a 
point. We do not yet have a topological interpretation of this torsion invariant, but hope to 
return to this soon. 

The main step in the proof of Theorem 11.11 is the construction of the heat kernel for the 
Hodge Laplacian on (M,g) using the methods of geometric microlocal analysis. This was 
carried out some time ago by Mooers [24] for spaces with isolated conic singularities, and the 
construction here is very similar in spirit (and most details), but must be done with careful 
attention to the possible terms which can arise in the asymptotic expansion of the pointwise 
trace. The 'geometric microlocal' method referred to here is the one pioneered by Melrose, 
see [22j, and developed by him and many others in the past few decades. It involves, for this 
problem, the careful study of the Schwartz kernel of the heat operator, as a polyhomogeneous 
distribution on a certain resolution of M x M x M + obtained by a sequence of real blowups. 
Heat trace asymptotics for simple edge spaces have also been studied by Briining and Seeley 
[8]; they approach these via the resolvent expansion, which although roughly equivalent is 
perhaps slightly less well adapted to the present purpose. 
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Our main theorem is a consequence of the following three main technical results. All 
notation is as above, but we also use some terminology which will be explained later. 

Theorem 1.2. For each degree k, the (rescaled) heat kernel H k '■= e~ l ^ k lifts to a blown up 
'heat space ' M% as a polyhomogeneous distribution. Its asymptotic expansions at all boundary 
faces of this space are determined by the indicial roots of Ajt- In particular, it is p~ m times 
a smooth function near the face td and p^^ 1 times a smooth function near the face ff ; where 
in either case p is a boundary defining function for the corresponding face. 

Theorem 1.3. The trace of the heat kernel H k is a polyhomogeneous distribution on R + with 
asymptotic expansion 

oo oo 

Tr H k (t) ~ jr +J2 C ^+E G ^ lo § f 

e=o e=o ee3 

as t — > 0, where 

3 = {leN \l + m-b even}. 

is the index set in the last sum. If g is an even asymptotically admissible metric, then this 
takes the simpler form 

oo oo 

Tr H k (t) ~ A et e ~ f + £ + 2 G ^ log '» 

e=o e=o eey 

where 3' = if (m — b) is odd and 3' = No = {0, 1, 2, . . .} if (m — b) is even. 

Theorem 1.4. Suppose that g^ is a family of asymptotically admissible metric which satisfy 
either conditions i) or ii) above. Then 

— \\.\\ RS -n 
H(m, 5m) -u. 



Our main result, the invariance of the analytic torsion under all deformations amongst even 
asymptotically admissible edge metrics, in case of the total and the edge dimension being both 
odd, has an important consequence towards the Ray-Singer conjecture for a certain class of 
stratified spaces. 

More precisely, consider a smooth closed Riemannian manifold (M m ,g) with a codimension 
two submanifold B m ~ 2 , and assume that m (and hence b = m — 2) is odd. View M as a simple 
edge space, and introduce a family of even asymptotically admissible edge metrics g a , where 
a £ (0, oo) is the cone angle in the direction normal to B, with g2 n = g. Since the Friedrichs 
extension of the Laplacian for (M, g^vr) coincides with the standard self-adjoint realization 
of the Laplacian on the smooth compact manifold (M, g) , Theorem 11.11 proves that for any 
a > 0, 

II \\RS _ || \\RS 
II • H(M, 5q ) - II • ll(M,g)- 

By the Cheeger-Miiller theorem, the right side is equal to the combinatorial Reidemeister 
torsion norm: 

|| \\RS _ || ||Reid 
II • \\(M,g) ~ II • llM • 

This gives a combinatorial meaning to the analytic torsion for the edge metrics g a . We 
might hope that the right hand side of this last equality agrees with the Reidemeister torsion 
norm defined in terms of the intersection homology of the stratified space (M,B), see |13| . 
Unfortunately, this last point is not so clear, see |12j . 
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This paper is organized as follows. We begin in £j2]with a more careful definition of the class 
of simple edge spaces, the various classes of incomplete edge metrics mentioned above, and 
an examination of the structure of the Hodge Laplacian for these metrics. We also consider 
the asymptotics of solutions to Aw = where w is in the maximal extension of this operator, 
and use this to characterize the Friedrichs extension. ^3] contains the parametrix construction 
for the heat kernel as an element of the calculus of heat operators on M\ we also introduce 
the even subcalculus and show that it contains when g is even. These results lead directly 
to a description of the asymptotics of the heat trace in §4; the proof that the analytic torsion 
is well-defined and the computation of its variation is given in $5j The appendix contains the 
proof of the composition formula for the heat calculus, as well as the even subcalculus. 
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2. Simple edge spaces and the Hodge Laplacian 

We begin by introducing the class of stratified spaces with simple edges and incomplete 
edge metrics, and then describe the structure of the Hodge Laplace operator on such spaces 
near the edge. We then define the class of polyhomogeneous distributions, and finally describe 
the Friedrichs extension of the Hodge Laplacian. 

2.1. Simple edge spaces. Let M be a compact stratified pseudomanifold, with top- 
dimensional (open, dense) stratum M, an open m-manifold, and a single lower-dimensional 
stratum B. There are a set of axioms, reviewed in [2], which regulate the differential topo- 
logical structure of a stratified space. For our purposes, there are two main consequences of 
these axioms. First, in our setting there are no higher depth strata on the frontier of B, so 
it is a compact smooth manifold; we could easily treat the case where B has components of 
varying dimension, but for simplicity we assume that B is connected, with dimB = b. Second, 
there is a neighbourhood IA of B in M, a 'radial function x defined onMnW and a smooth 
projection (j) : U — > B, which is a submersion on U n M, such that the preimages 4>~ 1 (q) 
are all diffeomorphic to truncated open cones C(F) over a compact smooth manifold F, with 
dimF = /, and where the restriction of the function x to each fibre 4>~ 1 (q) is a radial function 
on that cone. We suppose that the level set {x = 1} corresponds to the 'outer' boundary of 
the neighbourhood U, and we denote this manifold by Y . Thus Y is smooth and compact, 
and is the total space of a fibration : Y —> B with fibre F. Figure 1 illustrates this cone 
bundle structure. 

The class of metrics g we consider on such a space are the ones which restrict, on each 
conical fibre, to be asymptotically conic. To be more precise, let us say that g is rigid if it is 
arbitrary away from the open set U, but that in ti it has the following form: there is a smooth 
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Riemannian metric f) on B and a symmetric 2-tensor n on Y which restricts to a metric on 
each fibre F such that 

g\ u = dx 2 + cj)*\) + x 2 k. 




Figure 1. The tubular neighbourhood U in M. 

Hence in U, the induced metric on each fibre 4>~ l {q) is an exact warped-product conic 
metric. In fact, it is necessary for us to work with a slightly more restricted class of metrics: a 
metric g is called admissible if it is rigid and in addition eft : (Y, g\ Y ) — > (B, f)) is a Riemannian 
submersion. Recall that this means the following: if p E Y, then T p Y splits into vertical and 
horizontal subspaces, TpY @T^Y , where by definition TpY is the tangent space to the fibre 
of cj) through p and TpY is the orthogonal complement of this subspace. The new condition 
is that the restriction of the tensor n to TpY vanish. 

More generally, the metric g is asymptotically admissible if g = go + h where go is admissible 
and \h\ go — > at B. We shall typically assume that h has a poly homogeneous expansion in 
powers of x (see §2.4), or even more strongly that it is smooth to x = 0. In order to describe 
this accurately, it is convenient to pass from the singular space M to its resolution M, the 
manifold with boundary obtained by 'blowing up' the stratum B. In concrete terms, this 
corresponds to replacing each conic fibre C(F) = [0, l) x x Fj ~ (where (0,q) ~ (0, q') for 
any points q, q' £ F) with the cylinder [0, l) x x F, and correspondingly, replacing the cone 
bundle neighbourhood U with the associated bundle of cylinders. The C°° structure on M is 
the natural one induced from the cylindrical fibres. Notice that dM = Y. This resolution 
process is described in greater detail in [21], see also [2]. 

The advantage of considering objects on M rather than M is that there is now a clear 
meaning for a function or tensor to be smooth to x = 0. Thus we now consider metrics 

g\ u = dx 2 + <fi*t)(x) + x 2 k(x). 

where h(x) and k(x) depend smoothly on x E [0, 1). We can always decompose any such g as 
go + h where go = dx 2 + 0*f)(O) + x 2 k(0). Then g is called strongly asymptotically admissible if 
it is smooth and if \h\ go = 0(x 1+b ). It is called even if f)(x) and k(x) are smooth as functions 
of x 2 rather than just x, and if in addition h has a smooth expansion with only even powers of 
x. As in the introduction, we call g strongly asymptotically admissible or even in case these 
conditions are satisfied, respectively. 

2.2. Edge operators. The paper |21j discusses in great detail the theory of elliptic operators 
on simple edge spaces. We review some of the relevant material in this theory. 

Consider local coordinates (x, y, z) on M near the boundary face Y, where x is the radial 
coordinate, y is the lift of a local coordinate system on B and z restricts to coordinates on 
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each fibre F. Now define the class of edge vector fields V e on M: these are the vector fields 
on this space which are smooth even at x = and which are tangent to the fibres of Y at this 
boundary face. In this local coordinate system, any element of V e can be written as a sum of 
smooth multiples of the basic generators xd x , xd Vi and d Zj , which we write as 

V e = Span C oo {xd x , xd yi , xd Vb , d Zl , d Zf }. 

These are the basic objects of the edge theory. Notice that if g is a smooth asymptotically 
admissible metric on M, and if V, W G V e , then x~ 2 g(V,W) 6 C°°; in particular, the length 
of any edge vector field decays at least like x as x — > 0. 

We next introduce the class of differential edge operators Diff*(M). By definition, L € 
Diff^(M) if it can be written locally as a sum of products of elements of V e , with coefficients 
in C°°(M). Thus 

L= Y a j,a^(.x,y,z)(xd x y(xdy) a d^, 

j+\a\+\P\<m 

with each a,j a p smooth up to x = 0. More generally, if L acts between sections of two vector 
bundles, then L has this form with respect to suitable local trivializations, where each aj <at p 
is matrix- valued. The operator is called edge elliptic if its edge symbol 

e a m (L)(x,y,z;C,V,0 ■= Y a j,aA x ^y^ z )^V a C 13 

j+\a\ + \/3\=m 

is nonvanishing (or invertible, if matrix- valued) , for (£, rj, Q ^ (0,0,0). This has an invariant 
meaning as a function on the so-called edge cotangent bundle e T*M which is homogeneous 
of degree m on the fibres. 

There is an entire zoology of objects associated to edge geometry, but in the interests of 
space, we refer to the papers cited above for more on all of this. 

The fundamental tool in the analysis of elliptic edge operators is the space of pseudodif- 
ferential edge operators **(M). These operators are standard classical pseudodifferential 
operators in the interior, but are adapted to the degeneracy structure of elements of Diff* and 
form a suitably broad class of operators so that, in many cases, an elliptic differential edge 
operator L has a pseudodifferential edge parametrix G such that both GL — I and LG — I 
are compact on certain natural function spaces. This is the main content of [21j. 

One key idea in this theory is the use of two different model operators, which provide 
fundamental information about the edge elliptic operator L in the parametrix construction. 
The first of these is the indicial operator, 

I(L)= a jAI3 (0,y ,z)(sd s yd^ 

j+\/3\<m 

acting on functions on M + x F, where yo is some fixed point on B; by taking the Mellin 
transform in s, this reduces to the indicial family 

J C (L)= Yl w(0,yo,z)C J d^ z , 

j+\/S\<m 

which is a holomorphic family of unbounded Fredholm operators on L 2 (F). Values of £ for 
which Iq{L) is not invertible are called indicial roots of L; the indicial family is a generalization 
of the resolvent (of, say, Io(L)), and the indicial roots then play the role of the eigenvalues 
of this operator. Using the ellipticity of L, it may be shown that the set of indicial roots is 
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discrete in C and that the corresponding solutions of Iq(L)<P(z) lie in C°°(F). The second 
model is the normal operator, 

N(L) = Yl Ha,p(0>Vo,z)(sdsY{sd u ) a dg. 

j+\a\+\P\<m 

Here (s, u) E M + x M. b are linear variables on a half-space, which should be thought of as the 
inward pointing normal space to the fibre of Y through (0, yo,z). This seems to have almost 
the same complexity as L, but since it is translation invariant in u and dilation invariant in 
(s, u) jointly, it can be reduced by Fourier transform and rescaling to an operator on R + x F 
which is only slightly more complicated than I(L) to analyze. 

One main result in the theory is that if N(L) is invertible (on some fixed weighted L 2 
space) for each yo £ B, then L itself is Fredholm on the corresponding weighted L 2 space 
on M. The indicial roots of L determine a discrete set of weights such that L does not have 
closed range when acting on the corresponding weighted L 2 space. In any case, the point is 
simply that the inverse of N(L) is the main ingredient in the construction of a parametrix for 
L. Something similar is true for the construction of a heat kernel parametrix for L. 

We conclude this general discussion by observing that although the Laplacian and other 
natural elliptic operators on a simple edge space with incomplete edge metric are not quite 
edge operators, they are very closely related and can be studied by this edge theory. More 
specifically, if g is in one of the classes of incomplete edge metrics above and A g its scalar 
Laplacian, then A g = x~ 2 L where L E DifL? (M). The analogous statement is true for the 
Hodge Laplacian provided we trivialize the form bundles appropriately. We have already 
alluded to the edge cotangent bundle e T*M. In the local coordinates (x, y, z), this bundle has 
a local smooth basis of sections 

dx dy*_ dy b x f 
, , . . . , , uz , . . . , ctz 

XX X 

The edge /c-form bundle e A k (M) is simply the k th exterior power of this bundle, hence is 
generated locally by /c-fold wedge products of these sections. The correct assertion, then, 
is that A flj fe, as an operator acting on sections of e A k (M), has the form x~ 2 L where L E 
Diff2(M; e A fc (M)). 

We omit the proofs of these facts; they can either be checked by direct computation or else 
inferred by 'naturality' of edge structures, cf. [2"T] . 

2.3. The Hodge Laplacian. We now consider the structure of the Hodge Laplacian of an 
incomplete edge metric. More specifically, we describe some part of the structure of the normal 
operator N(x 2 A p ) and of its indicial roots. The material here is drawn from [T7], which in 
turn summarizes and presents in unified fashion various results proved by Bismut-Cheeger 
and Dai. 

The normal operator for x 2 A p at any yo E B acts on p-forms on the model edge M. b x C (F) = 
F with incomplete edge metric g = ds 2 + s 2 k(0) + |dn| 2 , and is naturally identified 
with s 2 times the Hodge Laplacian for that model metric. 

The first step is to consider the structure of the operator induced on the hypersurface 
S = {s = l}=I ll xF. For this, note that T(R b x F) splits into the sum of a 'vertical' and 
'horizontal' subspace, where the first is the tangent space to the F factor and the second is the 
tangent space to the Euclidean factor. This splitting is orthogonal, and induces a bigrading 

e A p (S)= e A J (lR 6 ) (g> A i (F) := e A^ e (S). 
j+e=p j+t=p 
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Let e Cl^ ,i (S) denote the space of sections of the corresponding summand in this bundle de- 
composition. The differentials and codifferentials on these factors satisfy 

d Rb : e n j ' e (s) -> e n j+1 '\s), d F ■. e n j ' e (s) e n j > t+1 (S) 

Sgb : 3 n j ' e {S) -> ^-^(S), 5 F : e n j ' e (S) -> e O^- 1 (5'). 

Furthermore, = d R b + dj? and £5 = (5 K 6 + 8p. It is not hard to check that the differential 
and codifferential induced on each level set S a = {s = a} has the form 

ds a = d R b + a~ l d F , 5 Sa = 5 Rb + a -1 ^. 

This can all be assembled into an expression for N(x 2 A p ) on R 6 xC(F). In order to simplify 
various calculations below, we write this operator using a rescaling of the form bundles, 
employed also by Br lining- Seeley [6] and in slightly different form in [17] . Thus, for each j, £ 
with j + £ = p, define 

<f> jjt : C°°(R+ W-*" 1 ^) e n j ' e (S)) -> e fF(M fe x C(F)), 

( r? ^)^/-l-//2 r?Ads + /-//2 /U) 

and denote by $ p the sum of these maps over all j + £ = p. It is not hard to check that 
(2.1) 

<3?p : L 2 J R + , L 2 ( e n^- 1 (S)® e n^(S),K(0) + \du\ 2 ),ds\ — ► L 2 ( e JF(M fe x C(F)),g) 
\ j+e= P J 

is an isometry, and a calculation yields 

(2-2) fc" 1 [s- 2 N(x 2 A p )] <S> P = (-^ + ^(A- 1/4)) + A B , 

where A is the nonnegative self-adjoint operator on A i ~ 1 (F) © A^(F) given by 

A*_i,f + (* - (/ + 3)/2) 2 2(-l)^, F 



(2.3) A 



2(-1)^_if A* F + (£-(/ + l)/2) 



If (M,g) is a simple edge space with an asymptotically admissible edge metric, then we can 
define a similar rescaling $ using powers of the defining function x in the neighbourhood U of 
B. Assuming that x is smooth onM\5 and equals 1 away from U, then this rescaling can be 
extended trivially to the rest of M. Conjugating by $ as before, then this rescaled operator 
Ap is a perturbation of (|2.2[) with higher order correction terms determined by the curvature 
of the Riemannian submersion (j) :Y —> B and the second fundamental forms of the fibres F 
in Y. We write this rescaled operator simply as A p if there is no danger of confusion. 

One reason for this transformation is that the indicial roots of A* have a particularly simple 
form: writing the eigenvalues of A as v 2 , with corresponding eigenform cpj, the corresponding 
indicial root of A p are the roots of the quadratic equation 

(2-4) - 7 (7 " 1) + v) ~ I = O if = \ ± ^1-1 + 4^ = \±vj, 

where we make the convention that every Vj > 0. If uq = 0, then the corresponding solutions 
are x l / 2 (f)Q and x 1//2 (log x)4>o- Note that s^~ Uj 4>j ^ L 2 when Uj > 1. 
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2.4. Polyhomogeneity. Before continuing, we must introduce the spaces of conormal and 
poly homogeneous distributions on a manifold with corners W. These are generalizations of 
and replacements for C°° functions in this setting, and are needed simply because solutions of 
edge elliptic and parabolic equations tend to have this form. Briefly, a function is conormal 
if it has stable regularity with respect to differentiations by arbitrary smooth vector fields 
which are tangent to all boundaries and corners of W; it is poly homogeneous if it has an 
asymptotic expansion at all boundary faces, and a product type expansion at all corners, 
in terms of powers of the boundary defining functions for the hypersurface faces, with all 
coefficients depending smoothly on the tangential variables. The exponents which appear in 
these expansions may be arbitrary complex numbers, and we also allow nonnegative integer 
powers of the logs of the defining functions as factors. 

To state all of this more formally, let W be a compact manifold with corners, with all 
boundary faces embedded, and {Hi}f =1 the set of all hypersurface boundaries of W. As part 
of our definition of manifold with corners, we require that each Hi be embedded, so it has a 
boundary defining function pi, i.e. pi is a smooth nonnegative function on W which vanishes 
simply on Hi and is strictly positive on W \ Hi. For any multi-index A = (Ai, . . . , Atv) G R-^ 
we write p x = p Xl . . . p^f . Finally, denote by V&(W) the space of all smooth vector fields on 
W which are unconstrained in the interior but which lie tangent to all boundary faces. Thus, 
if p is a point on a codimension k corner of W, then there exists a local coordinate chart near 
p, (x±, . . . ,Xk,yi, ■ ■ ■ ,y n -k) with each Xj G [0, e) and yj G (— e, e), and very similarly to our 
definition of V e in §2.2, 

V b (W) = Span C oo {xid X] ,d ys , i, j = 1, . . . , k, s = 1, . . . , n - k}. 

Definition 2.1. A distribution w on W is said to be conormal of order A if it is of stable 
regularity with respect to elements ofVb(W), i.e. if it lies in the space 

A X {W) = {w G p x L°°(W) :V 1 ...V e we p x L°°{W) VVj G V b {W),£ > 0}. 

We also write A*{W) = (J\A X (W). Note that any conormal distribution is smooth in the 
interior ofW. Furthermore, L°° could be replaced by any other fixed space, e.g. some (poly- 
nomial^) weighted L 2 or LP space; the union of these spaces over all A is the same in any 
case, only the weight and regularity scale would change. 

Next, an index set E = {(j,p)} cCxNis the set of exponents associated to an expansion 
at the face Hi. We require that it satisfies the following hypotheses: 

(i) Each half-plane Re ( < C contains only finitely many 7; 

(ii) For each 7, there is a P{i) G No such that (j,p) G Ei for every < p < P("f) < 00; 

(iii) // ( 7 ,p) G Ei, then (7 + j,p) G Ei for all j G N. 

An index family £ = (E±, . . . , En) is an N -tuple of index sets, where Ej is associated to the 
face Hj ofW. 

Finally, we say that w is poly homogeneous on W with index family £ if it has expansions 
at the various boundary hyper surf aces with exponents determined by the index family £: 

A% hg (W) = {w G A*(W) : w~ ofyplOogPiT near Hi}, 

(7,p)eE< 

where all coefficients are themselves polyhomogeneous, G •A.^ g {Hi). The index family f W 
is the obvious one induced from £ at the boundary faces Hi n Hj of Hi. 



Rafe Mazzeo and Boris Vertman 



11 



It is a simple consequence of this definition that if w is poly homogeneous, then it has a 
product type expansion at any corner D ... Pi of W of the form 

w ~ ^a 7iP p 7 (logp) p , 

where now 7 and p are multi-indices, and with coefficient functions conormal on that corner. 
We recall also that by saying that the expansion for w is asymptotic, we mean that the 
difference between w and any finite portion of the expansion vanishes at the rate of the next 
term in the expansion, with a corresponding property for all higher derivatives. 

In the next subsection we shall also encounter distributions which are conormal and 'par- 
tially polyhomogeneous', i.e. they are conormal and have a finite expansion up to some order 
of decay, with a remainder term which is only conormal. We do not introduce special notation 
for these spaces. 

We refer to [21] for more more details about poly homogeneity. 

2.5. The Friedrichs extension. Since the compact simple edge space (M,g) is incomplete, 
the Hodge Laplacian may not be essentially self-adjoint on the core domain C^ D £l p (M), so 
we must consider how to impose boundary conditions at the edge to obtain closed, or even 
better, self-adjoint, extensions. For spaces with isolated conic singularities, this was first 
accomplished by Cheeger [10]. Further and more systematic studies in the conic setting for 
general Dirac-type operators appear in [20] and [15], and see [24], [E] for results about the 
associated heat equation. That setting is tractable because the extension problem is finite 
dimensional. When the edge has positive dimension, the requisite analysis is more intricate. 
We restrict attention to the Friedrichs extension since this requires much less machinery to 
define and characterize. 

Let A* = A p denote the rescaled Hodge Laplace operator acting on differential forms of 
degree p on the compact simple edge space M. Consider the space of L 2 forms L 2 Cl p (M), 
with respect to any choice of (polyhomogeneous) incomplete iterated edge metric on M, as 
well as the associated edge Sobolev spaces 

H e e tt p (M) := {u G L 2 QP(M) \ V\ ■ ■ ■ VjU G L 2 Q P (M) for Vi G V e and for any j < £}. 

We often use to denote the intersection over all t of the spaces Hg . The maximal domain 
of A p is, by definition, 

£>max(A p ) := {u G L 2 n p (M) I A p u G L 2 n p (M)}, 

where A p u G L 2 is initially understood in the distributional sense. Similarly, the minimal 
domain of A p is defined as 

£> min (A p ) := {u G £> max (A p ) | 3uj G C£°JF such that 

Uj —7- u and A p uj — > A p u both in L 2 Q P }. 

These are the domains of the maximal and minimal extensions of A p on the core domain 
C^°rP(M). The set of all closed extensions of A p is in bijective correspondence with the 
closed subspaces of the quotient T> max /V m j n ; furthermore, since A p is symmetric on the core 
domain, self-adjoint extensions are in bijective correspondence with the subspaces of this 
quotient which are Lagrangian with respect to a certain natural symplectic form induced 
from the boundary contributions in an integration by parts formula, see pOj and [15J. 

This motivates the problem of characterizing elements in the maximal domain. For spaces 
with isolated conic singularities, this is straightforward and leads to an explicit parametriza- 
tion of all closed and self-adjoint extensions of A p . For spaces with simple edges singularities, 
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however, the complete characterization remains functional analytic, although there are many 
explicit extensions which parallel the definitions in the conic setting. Different choices of 
closed extensions in either case correspond to what are sometimes called ideal boundary con- 
ditions at the singular stratum. We begin with a result which holds in both the conic and 
edge settings. 

Lemma 2.2. Let (M,g) be compact with simple edge singularity. Then 
^max(A p ) C H™n p (M)+x 2 H 2 e n p (M) C H 2 n p (M), 
^min(Ap) C x 2 ~ € H 2 Q p (M) for any e > 0. 

If M has only conic singularities, then any w G T> max (A p ) admits an asymptotic expansion 

A i (x^-^aT(z), Vj^O 

(2.5) w ~ w + J2 x * +Vja j( z ) + { i w£x 2 H 2 n p , 

where the numbers l/2±Vj are the indicial roots ofA p with Vj G [0, 1) and each a^(z) G C°°(F) 
is a solution of the corresponding indicial operator. If M has a simple edge singularity, then 
w admits an analogous expansion with coefficients a^(y,z) now depending on both y and z, 
but this is an asymptotic expansion only in a weak sense, i.e. there is an expansion of the 
pairing J B w(x,y, z)x(y) dy for any test functions on B. 

Remark 2.3. We refer to |21| for a careful explanation of such weak expansions. It is the 
failure of the expansion to hold in a strong sense which makes this result more difficult to use 
for simple edge spaces with dim B > 0. See also |17| and [2] for the case where there are no 
indicial roots in the critical range, so that £> max = ^min- 

Proof. Consider any w G L 2 Q P (M) with A p w = f G L 2 Q P (M). Then x 2 A p w = x 2 f G 
x 2 L 2 Q p (M), which is useful since the operator x 2 A p is an elliptic differential edge operator 
of order 2, which we denote by C. Applying [21L Theorem 3.8], we obtain a parametrix 
B G v & f T 2 (M) for C in the small edge calculus, which satisfies 

BC = I-K, K£ <&-°°{M). 

Corollary 3.23 in [21] asserts 

B : x s H ( n p (M) -)• x s H e+2 n p (M), 

n ■. x s H e n p (M) x s H™n p {M), 

for any s G M, £ G N. Applying B to Cw = x 2 f yields w = B(x 2 f) + 1Z(w), and this lies in 
H%° + x 2 H 2 as claimed. 

If M has only isolated conic singularities, then there is a considerable sharpening of the 
parametrix construction, see |21| or [22]. It is then possible to choose a much better parametrix 
B which has the property that the remainder term TZ = I — BC maps an arbitrary element 
of L 2 Q P (M) into a polyhomogeneous distribution. We do not describe that construction, but 
it is very similar in spirit (and simpler than) the construction for the heat parametrix in §4 
below. In any case, granting this, then the corresponding equation u = B{x 2 f)+7tw now gives 
that u G x 2 H 2 {M) + *4 p h g , i.e. it has a partial expansion up to order x 2 as in the statement 
of the lemma. Once we know that u has such a partial expansion, we can determine the 
structure of the terms which vanish less quickly than x 2 simply by substituting this expansion 
into the equation and calculating formally. The only possible terms x^^z), 7 < 2, for which 
C(x' y (j)(z)) = 0(x 2 ) are those with 7 indicial and <f> a corresponding solution. 
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The proof that such a partial expansion holds in the weak sense in the simple edge case 
uses the Mellin transform, and is explained carefully in §7 of |21| . 

If w G T> max , then a straightforward calculation shows that any function which can be 
approximated in the graph norm cannot have any terms of the form x 7 a(y, z) (even if just in 
a weak expansion). This proves the claim about the minimal domain. □ 

To characterize the Priedrichs extension of the Hodge Laplacian, observe that A = Q) p A p 
factors as A = DoD, where D = d+5 is the Gaufi-Bonnet operator on Vt*(M). Letting D act 
on sections of e Q*(M), as above, then xD is an elliptic edge operator of order 1. Its normal 
operator N(xD) yo at any yo G Y acts on the model edge R b x C(F) and in fact s~ l N(xD) yo is 
naturally identified with the Gaufi-Bonnet operator Dg for the model incomplete edge metric 
g = ds 2 + s 2 k + \du\ 2 . Conjugating by the unitary transformation $ from (12. ip we calculate 
that, according to the splitting fT = fi cvcn fi odd , 



(2.6) D 



g ~ \f 8 , l 



ds 



(-& + $P)+D* 

°)+D Rb , 

+ £ (A+ - 1/4) + A Rb 
" I _^ + ^ ( ^__;i /!) + a, 

6 



and so 

/_A 2 i i 
(2.7) Dl= 9s 



Here D K b is the Gaufi-Bonnet operator on R , P is a self-adjoint first order differential operator 
on n*(F), and 

(P + l/2) 2 = A+, 
(2.8) 9 
V ; (P-1/2) 2 = A_. 

If (M,g) is a simple edge space with an admissible edge metric, then D is a perturbation 
of (d s + s~ 1 P) + Dy with higher order correction terms determined by the curvature of the 
Riemannian submersion <j) : Y — )• B and the second fundamental forms of the fibres F in 7. 

Just as for the Hodge Laplacian, we can define the minimal and maximal extensions, D ni \ n 
and -D max of D; their domains are V(D m \ n ) = V m \ n (D) and V(D ma _ x ) = V(D max ). We now 
characterize elements in these domains using arguments similar to those in Lemma [ 



Lemma 2.4. Let (M, g) be compact with simple edge singularity. Then, for any e > 

^max(-D) C H™n*(M) +x 1 - e HlQ p (M) C HlQ*(M). 

Any w £ I'max(-D) admits a weak asymptotic expansion 
N 

w ~ J2 x ~ Uj+1/2b j(y> z ) vj € Spec(A+®A-) n (0,1) w G x 1 ~ e H^S}* . 

T/ie coefficients bj(y,z) are smooth in z, but typically have negative Sobolev regularity in y. 
On the other hand, 

T> min (D) = {weV max (D) | bj(y,z) = 0,j = 1,...,N}. 

Proof. Consider w = (w+,W-) G T> max (D) with respect to the decomposition of D into even 
and odd components. The existence and nature of the terms in the weak asymptotic expansion 
for elements w± is derived exactly as for the Laplacian in Lemma 12.21 see |21} Theorem 7.3]. 

The exponents 7 in the expansion of w+ arise as indicial roots of (sd s + P), in the 
interval (—1/2,1/2). Consequently 7 G — Spec (P) D (—1/2,1/2), with the coefficient 
&+(— 7) G T)'(Y, E(— 7)) being a distribution in Y with values in the finite dimensional 
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(— 7)-eigenspace of P. Since on forms of even degree, (P + 1/2) 2 = A+, we see that 
7 = —v + 1/2, v 2 G SpecA + n (0, 1). Hence the expansion of w + takes the form 

w + ~ x76 +(-7) + w+, w+ G x 1_e .H£fi*. 

Similarly, the exponents 7 in the expansion of u>_ arise as indicial roots of (— sd s +P), in the 
interval (-1/2,1/2). Consequently 7 G Spec (P) n (-1/2,1/2), with the coefficient 6.(7) G 
T)'(Y, E("f)) being a distribution in Y with values in the finite dimensional 7-eigenspace of 
P. Since on the odd components (P — 1/2) 2 = yl_, we find again 7 = —u + 1/2, v 2 G 
SpecA_ D (0, 1). Hence the expansion of W- has the form 

w - ^^x^b^i^ + w-, w- G x^Hln*. 

It remains to establish the claim on -D mm . Note that if /3±(7) G C°°(Y, E(j)) with 7 G 
Spec(P)n (-1/2,1/2), then (x^P + (j), x^P-d) G V mSuX (D). Now consider to = (w + ,W-) G 
^min(-D) C T> max (D) with asymptotic expansion as above, and suppose that v G T> max (D) is 
given by a sum of test sections (x _7 /3+(7), x 7 /3_(7)), cut off smoothly at x = 1 and extended 
to all of M. Then 

(2.9) (Dw,v) - (w,Dv) = £(& + ( 7 ),0_( 7 )) + (6_( 7 ),/3 + ( 7 )), 

7 

where the summation is over 7 G Spec (P) n (—1/2,1/2). This vanishes because and v 
lie in adjoint domains, namely V m \ n and 2? m axj respectively. Since the coefficients /3±( 7 ) are 
arbitrary, we deduce that each 6-j- (7) = as a distribution. Hence the vanishing of these 
conditions is at least a necessary criterion to lie in the minimal domain. It is also sufficient 
since any v G T>(D max ) can be approximated in L 2 by polyhomogeneous forms, and if all 
coefficients for w vanish, then {Dw,v) = {w,Dv), and this is sufficient to show that w lies in 
the domain adjoint to -D m ax> i.e. in V(D m \ n ). □ 

Finally, recall the characterization of the Friedrichs extension for the Hodge Laplacian 

A Fr = AnaxAnin, T>(/SF T ) = {w G V min (D) : Dw G AnaxP)}- 

Thus, using the notation of Lemma 12.21 Lemma 12.41 implies the following. 

Proposition 2.5. Let (M,g) be compact with simple edge singularity and an adapted incom- 
plete edge metric g, and the Friedrichs extension for the Hodge Laplacian. Then 

P(A Fr ) = {w G IWA) I aj(w;z)=0,j = 1, ...,#}. 

This result is well-known especially in the conic case. Henceforth we omit the superscript 
Fr since we deal exclusively with the Friedrichs extension. 

3. Construction of the heat kernel 

In this section we develop the 'heat calculus' on a simple edge space M and prove that 
the heat kernel of the Friedrichs extension of the Hodge Laplacian lies in this calculus. This 
construction follows a now standard path, and in fact is almost identical to the construction 
in [21] of the heat kernel for the Laplacian on spaces with isolated conic singularities; other 
related constructions appear in [22], [I], and see also [3] for an expository account of several 
other heat kernel constructions using these same techniques of geometric microlocal analysis. 

The heat kernel is a priori a distribution on [0, 00) x (M) 2 which is C°° on (0, 00) x (int M) . 
The general idea is to study its singular structure by lifting it to a 'heat space' M^, which is a 
resolution of M + x (M) 2 obtained by blowing up certain submanifolds of the boundary of this 
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space; elements of the heat calculus are, by definition, operators with Schwartz kernels which 
are polyhomogeneous on M? . An iterative parametrix construction is used to construct a good 
approximation to the fundamental solution operator for dt + A g in this class of operators, and 
a standard regularity argument shows that the true heat kernel itself lies in this calculus too. 
Careful bookkeeping yields very precise information about the small-time asymptotics of this 
kernel near the singular strata, which leads eventually to the main result of this paper. 

3.1. The heat double space. Fix an adapted coordinate chart (x, y, z) on M near the edge, 
or equivalently, on M near its boundary. Taking two copies of this chart as well as the time 
coordinate t yields a coordinate system (t, x, y, z, x, y, z) on M + x (M) 2 , valid near the diagonal 
(x,y,z) = (x,y,z). We define the resolution in two steps. The first is to blow up the fibre 
diagonal of dM at t = 0, i.e. the submanifold 

FD = {(0,0,y,z,0,y,5) :y = y}, 

using the parabolic homogeneity of the problem. The space which results from this is denoted 
FD, {(it}]; this notation indicates both the submanifold that is to be blown up (FD) and 
the direction dt which is scaled parabolically. We refer to [13], see also [3J, for a more careful 
discussion of parabolic blowups. This space is pictured in the following figure. 




tf 



Figure 2. The first blowup of the heat space. 

It has four codimension one boundary faces, which we denote ff, If, rf and tf (for 'front 
face', 'left face', 'right face' and 'temporal face'), respectively. One way of understanding this 
blowup is that it is the smallest manifold with corners on which the lifts of smooth functions 
from M + x (M ) 2 and the 'parabolic polar coordinates' R = \fx 2 + x 2 + t } = (t/p 2 ,x/p,x/p) 
are all C°°. These polar coordinates are difficult to use in computations, so we typically work 
with projective coordinate charts instead. Thus away from tf we use 

r x ~ x y — y ~ ~ 

(3.1) P = Vt, £ = £ = w = — /=-> z > 2/> z - 

Here p, £ and £ are defining functions for the faces ff , rf and If respectively, away from tf . On 
the other hand, away from If we use the chart 

t x y — y _ _ _ 

(3.2) r = ^r, s = -, to = , z, x, y, z, 

X z X X 

so that t,s,x are defining functions of tf , rf and ff , respectively. There are also analogous 
coordinates valid away from rf , obtained by interchanging the roles of x and x. 

Now, let Dq = {(0, x, y, z, x, y, z)} be the diagonal of (M) 2 at t = and D\ its lift to this 
intermediate heat space [M + x (M) 2 ; FD, {dt}]. In the two coordinate systems (|3.ip and (|3.2p 
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above, 

Di = {p = 0,^ = ^,w = 0,z = z} = {t = 0,s = 1,uj = 0,z = z}, 

respectively. The second and final step is to blow up [M + x (M) 2 ; FD, {dt}] along D\, again 
parabolically in the t direction. This yields the final heat space: 



Ml 



[R+ x {M) 2 ;FD,{dt}];D 1 ,{dt} 



This has the same four boundary hypersurfaces as before, but also a new one, denoted td (for 
temporal diagonal), created in the final blowup. This is illustrated in the picture below. 




Figure 3. The blown up heat-space M%. 

Using two pairs of interior coordinates (x, y, z) and (x, y, z) on M, then projective coordi- 
nates near td and away from ff are given by 

(3.3) T = Vt, S = ((x - x)/T, (y - y)/T, (z - z)/T), x,y,z. 

For a coordinate chart near td valid up to the corner td D ff but away from tf and If, use the 
previous projective coordinate system to define 

s — 1 z — Z IS) 

(3.4) 77 = a/t, a = , C = , M = — , x > V, z - 

rj r] T) 

Thus T and 77 are the defining functions for td in these latter two coordinate systems, while 
a; is a defining function for ff in the first of these projective coordinate systems. Away from 
td, it is sufficient to use the projective coordinate systems from the intermediate space. 
There is a canonical blow-down map 

/3 : M% — >R + xM 2 

which can be written explicitly in any of the local coordinate systems above. It is pertinent 
that the restriction of (3 to td is a fibration, with each fibre /? _1 ((0, w , w)) a closed 'parabolic' 
hemisphere and base the lifted diagonal of (M) 2 , which is diffeomorphic to a copy of M. 
Similarly, the restriction of (3 to ff is a fibration over the diagonal of dM with fibre equal to 
the product of a 'parabolic' quarter-sphere and the fibre of M, i.e. S*™^ 1 x F. 

3.2. Heat calculus and parametrix construction. We now define the heat calculus on 
M. This is a set of operators, characterized through the polyhomogeneity properties of the 
lifts of their Schwartz kernels to the double heat space M 2 . Recall that a function Ka on 
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M + x M 2 acts as an operator in one of two ways: first, it carries functions on M to functions 
on R+ x M by 

0i — ► lK A (t,z,z)(j)(z)dz, 
JM 

and second, it acts on functions on R + x M by 

/i — > lK A (t-s,z,z)f(s,z)dz. 
Jo JM 

Definition 3.1. Let £ = (E\f,E T f) be an index family for the two side faces of Ml. We define 
ty'\(M) to be the space of all operators A with Schwartz kernels Ka which are pushforwards 
from poly homogeneous functions Ka on M?, with index family {(—b — 3 + £ + j,Q) : j G N} at 
ff, {(— m + p + j,0) : j > 0} at td, at tf and £ for the two side faces of M%. When p = oo, 
Etf = 0- For simplicity we usually denote the lifted Schwartz kernel simply by Ka again. 

This is called a calculus because, at least under certain restrictions on the index sets at the 
side faces, one has that 

where £" is a new index family constructed explicitly from E and E' . We do not need the 
most general form of this composition rule, but only the special case when p = p' = oo, which 
is easier to prove. We defer the rather technical proof to the appendix. 

For simplicity of notation in this section, we define the heat calculus and carry out the 
parametrix construction only for scalar operators. The changes needed to handle operators 
acting between sections of a vector bundle V, in particular V = A P T*M, are minor; one must 
simply tensor at the appropriate places with V or End(V). We leave details of this to the 
reader, and shall reintroduce the bundle notation when it is needed more specifically. 

The main point of this section is that it is possible to construct a parametrix for the solution 
operator of C = dt + A as an element in VPg!?^ (M) for some index family £; indeed, using 
this parametrix one can then show that the exact solution operator itself lies in this same 
space too. We are only interested here in the case where A = A^ r , but the same construction 
applies equally well for the heat equation associated to many other natural geometric operators 
associated to an incomplete edge metric on M. This parametrix construction follows the 
scheme in [23], as explicated further in [3]. 

The construction is a standard one in geometric microlocal analysis: it proceeds by positing 
that the heat kernel lies in (M) for some £, and then solving for the successive terms 

in the expansion at certain faces of M% using the restrictions to the these faces of the lift of 
tC This leads to a parametrix for which the remainder vanishes to high, or even infinite, 
order at these faces. A better parametrix, which vanishes to infinite order at other faces, is 
constructed using the composition formula, and a short argument then yields that the true 
heat kernel is an element of this heat calculus and that the parametrix captures the form of 
the asymptotic expansions at all boundaries. 

It is helpful to work with the operator tC rather than just £, simply because it lifts to 
an operator which is smooth at the faces ff and td of M^; this change is unimportant since 
CH = implies tCH = 0, etc. To commence, then, consider the lift of tC from R + x M, 
where A acts on the left copy of M, to Ml. Thus locally near ff , td and tf, tC is a 6-operator 
on M|, i.e. a sum of products of elements of Vfe(M^). (However, near rf it can be put into this 
form only by multiplying by r 2 ). This is easily checked using the local projective coordinate 
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systems above. For example, near ff we compute that td x = pd^, td% = <9|, and so forth. 
Because of this structure, it makes sense to restrict tC to one of these boundary faces. We 
are particularly interested in its restrictions to td and ff , which we call the normal operators 
of tC at these faces, and denote by N t d(tC) and Ng(tJ0), respectively. Using the projective 
coordinate system (T, S,£) near td, we compute that 

(3.5) N td (t£) = l -Td T + A H - \ £ 3^ , 

while using the projective coordinates (r,s,w,z,x,y,z) near ff we see that 

(3.6) N(tC) =r(d T - d 2 s + s~ 2 (A - 1/4)) + A* 1 ) , 

where A is the rescaled operator (|2,3p acting in the variable z. Note that Ng{tC) acts tan- 
gentially to the fibres S™^ 1 x F of ff , i.e. it involve no derivatives with respect to (x,y,z). 

We now search for an element H 6 V? 2 ^ (M) for which tCoH vanishes (to as high order as 
possible) at as many faces of M 2 as can be arranged. The index family £ will be determined 
in the course of the construction. 

The expansion of H near td has a 'universal' form. Indeed, the expression for N t( ^(tC) above 
does not depend on the edge geometry, but is just the same as for the corresponding operator 
on a closed manifold without edges or other singularities. Using the projective coordinates 
T = y/i and 3 = — /T, write H ~ ^ T~ m+ ^ Hj near this face. The reason we have chosen 
the leading exponent to equal — m is to be compatible with the initial condition 

H\ t=0 = 6(t-0 = T-™6(Z). 

Applying N t d(t£) to this formal expansion for H shows that 

(A S -^3.a s + ^)^ = Q„ 

where Qj is an inhomogeneous term depending on Hq, . . . ,Hj_i. In particular, Qq = and 
H = (47r)- n / 2 T- n exp(-|3| 2 ). One can verify inductively that each Hj decays rapidly in 3, 
and moreover, that this construction is uniform in the parameter £ = (x,y,z), even up to the 
front face where x = 0. We refer to [22] and [3] for more leisurely treatments of this step. 

We can also solve away at least the leading order coefficient of H at ff. As in the previous 
case, expanding H near this face and writing Nff(H) for its leading coefficient there, we see 
that 

N s {tC oH) = N B (t£) o N B (H). 

To make this vanish, it is reasonable to choose Ng(H) to equal the fundamental solution for 
the heat operator Ng(tC). This heat operator is simply r times the heat operator for the 
Friedrichs extension of the Hodge Laplacian on the product C(F) x M. b . In other words, we 
set 

N S (H) = H c ^ (r, s, z, 1, z)H Rb (r, <j, 0). 
The first factor here is the Friedrichs heat kernel of the rescaled operator 

< 3 - 7 > a<w ~£ + ?H)' 

as in [9], |24| and |20j . Note that Ng(H) acts on the fibres of ff. It has a singularity at 
r = 0, s = 1, oj = in ff, which is precisely the intersection ff (ltd. This matches the singularity 
of Ntd(H) along td at this same intersection; this follows from the fact that N t ^(tC)HQ = 
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has a unique tempered solution up to scale, and we have chosen the normalizing constant to 
make the value at T = equal to the delta function. We refer to [3j for more on this point. 

We now find an element G C^O with this prescribed asymptotic data at these two 

faces. First choose some function H' near td which has the given sequence of functions Hj as 
the coefficients in its asymptotic expansion as T — > 0. Fix a smooth cutoff function x{ a ) which 
equals 1 for a near and which vanishes for a > c > and regard x(Pfr)% 1_b -^%(-£0 as being 
defined on M%. By the compatibility discussed above, we see that these two requirements for 
fit together at ff n td. 

This discussion provides most of the proof of the 

Proposition 3.2. There exists an element € SSf^\ (M), where Eif = E v { is an index set 
determined by the set of indicial roots of x 2 A®, i.e. the numbers | + Uj, E spec (A), such 

that tCHM = P (1) G ^l'^ £{1 \M), where = {E Yi ,E ri - 1), and with 
lim H^\t, x, y, z, x,y,z) = 5(x — x)5(y — y)5(z — 2). 

Proof. The only point which is not immediately apparent from the construction of H^ l > above 
is the characterization of its expansions at If and rf. These are inherited directly from the 
expansions of the Friedrichs heat kernel on the cone, and so we turn to a closer examination 
of this kernel. 

The Hodge Laplacian on the complete cone C(F) can be described using separation of 
variables. Following [20, Prop. 2.3.9], and using the eigendecomposition of the fibre operator 
A, Ap separates into the ordinary differential operators 

I ^-f- + l(u 2 -- 
dx 2 x 2 \ 4 

,2 



where v G spec (A). There is an explicit formula for the corresponding Friedrichs heat kernel: 

(3.8) e-^(x,x) = \ % ixxfl 2 I v ) e"^, 
and this lifts in projective coordinates to 

(3.9) H c ^(t, s, z,1,z) = x~ l ^V~s ■ h (£) e-^M^Mz), 

where <f> v is the eigenform associated to v 2 . Using standard estimates for the eigenvalues and 
eigenforms of A and the classical bounds for the functions /„, we can prove that this sum 
converges locally uniformly in C°°. We then write 

(3.10) H C(F)xR» g) 1? ^ Wj ~ Q) = H C(F) (T) 8 ^ Zf uz)H ^ (T) Wj q)) 
where 

F Rb (r,^,0) = x- b -^—- h e-^ 2 ^. 



(4vr)V2 



Since M + x C{F) x M. b is naturally identified with the interior of each fibre of ff, then directly 
following |24| we see that H c ( F ) x ^ b [ s poly homogeneous on the front face of M?. Furthermore, 
the index sets of H^ 1 ' at the left and right faces are exactly the same as those for H c( - F \ The 
index set at If, i.e. as s — > 0, can be read off directly from f|3.9|) : the exponents are simply the 
indicial roots v > 0. By symmetry of the heat kernel, the index set at rf is exactly the same. 
On the other hand, the index set for the remainder term at rf is lowered by one because 
the operator tC lowers terms by two orders but the leading term is killed. □ 
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The next step in this construction involves choosing a slightly finer parametrix which 
leaves an error which also vanishes to infinite order along the entire left face as well. This is 
done by an iterative construction similar to, but even simpler than, the procedure above at 
td. 

Proposition 3.3. There exists an element £ ^I'^(M), where £ is the same index 

family as above, such that t£ff( 2 ) = P^ £ $Z,°£ E if,°° (m) and lim^ # (2) = Id. 

Proof. This step proceeds exactly as in [24]. The error term P^ from the previous step has 
an expansion along rf, and we use tC to subtract off the successive terms of this expansion. 
Let J denote a kernel which has asymptotic sum at this face with terms equal to the ones 
obtained in this procedure; note that we can assume that J vanishes to first order at ff and 
to infinite order at tf. We see that = H^- 1 ' + J has the all the desired properties. 

In order to eliminate a term s 7 a in the asymptotic expansion of the error term P^ at rf , 
it is only necessary to solve the indicial equation 

(-d 2 s + s~ 2 {A - l/4))u = s 7 (r _1 a) . 

This is because all other terms in the expansion of tC at rf lower the exponent in s by at most 
one, while the indicial part lowers exponent by two. Note that r, u, x, y and z only enter this 
equation as parameters. We have already discussed how to solve this equation on the cone 
C(F) using the Mellin transform. The solution is polyhomogeneous in all variables, including 
the parameters. It has leading order 7 + 2 at rf . 

Iterating this argument, asymptotically summing these solutions, and adding the re- 
sulting kernel to we obtain a new parametrix 

H (2) £ grW^jif) with the 

same in- 
dex set £ = (E\f,E T f) at the right and left boundary faces, and with a new error term 

p{2) G ^oo,E l{ ,oo^ M y u 

We have now constructed a parametrix such that tCH^ = P^ vanishes to infinite 
order along If and all along the bottom faces tf and td, and which vanishes to first order at 
ff . We can regard any of these kernels as acting on functions on R + x M in the usual way, by 

(Ku)(t,x,y, z) = / / K(t — s,x,y,z,x,y,z^u(s,x,y,z^dsdxdyd'z. 
Jo Jm 

The identity operator corresponds to the kernel = 6(x — x)5(y — y)5(z — z). Viewed in 
this way, 

tCH^ =U + PW. 

The final stage in the parametrix construction is to consider the formal Neumann series 

00 

(id + p^))- 1 = id + Y,(- p(2) ) j ■= id + p(3) > 

i=i 

and then define 

=HW(Id + P&). 

In fact, a slightly finer analysis, see |22| and also [5], shows that this formal series is actually 
convergent. We refer to these sources for the necessary estimates, both in the case of compact 
manifolds and for manifolds with cylindrical ends, but it is easily seen that everything there 
transports immediately to this setting (and many other ones as well). Indeed, this is a general 
feature of such Volterra series. 
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The composition formula for the heat edge calculus, which we discuss in the appendix, 
yields the following: 

(p( 2 ))J £ ^3j,oo,B lf ,oo^ p(3) ^ ^oo,oo,£ lf ,oo^ fj(2) p(S) g ^00,00,-%,^ 

In other words, successive composition of P^ 2 ) with itself produces an operator which vanishes 
to higher and higher order on ff. The final operator satisfies tCH^ = Id in the sense 
of operators on M+ x M; equivalently, tCH^ = as an operator from functions on M to 
functions on R + x M, and # (3) = Id at t = 0. 

We have now produced a kernel i? = .ff( 3 ) which has the following three properties: 

• tCH = 0; 

• H\ t=0 = S(x - x)8{y - y)S(z - z); 

• The range of H lies in the Friedrichs domain of C for all t > 0. 

These three properties characterize the Friedrichs heat kernel uniquely. Hence H = . 
We have proved the 

Proposition 3.4. Let {M,g) be a compact simple edge space with asymptotically admissible 
metric g. Then the heat kernel for the Friedrichs extension of the (rescaled) Hodge Laplace 
operator on M is an element of the heat space ^\ (M;A k M), where 

Ex = E r{ = {u + i + £ : v 2 E spec (A), i G N }. 

3.3. The even subcalculus. The parametrix construction we have given already contains a 
significant amount of information about the asymptotic properties of the heat kernel at each 
of the boundaries of M?. We now explain how, if one makes stronger assumptions about the 
asymptotic structure of the metric g near the edge, one obtains stronger conclusions about the 
structure of these asymptotics. This will be important in the next section when we examine 
the heat trace and its asymptotics. 

To be more specific, we shall show that if the asymptotically admissible metric g is even, 
then the Friedrichs heat kernel lies in a distinguished subalgebra of the heat calculus which 
we call the even subcalculus. Elements of this even subcalculus are characterized by certain 
parity conditions on the coefficients in their asymptotic expansions, near both the faces td 
and ff of M?. As we exploit later, because of these parity conditions, many terms in the heat 
trace expansions of such elements vanish. Conditions of this type appear in many places in 
the literature, see [28] for example, for a discussion of the closely related Kontsevich-Vishik 
canonical trace on a similar subclass of pseudodifferential operators on a closed manifold. 
Even more closely related to the discussion here, Melrose defines an even subcalculus for heat 
operators on closed manifolds and on manifolds with asymptotically cylindrical ends [22], and 
we refer also to the even subcalculi in the papers by Albin [H Cor. 8.9] and Albin-Rochon [1] 
in other geometric settings. 

We shall impose separate 'evenness' conditions on the asymptotic expansions of elements 
of ^*_ h at the two faces td and ff . The condition at td is exactly the one introduced in [22] 
(and which appears at least implicitly in many other places), while the one at ff is directly 
inspired by this but tailored to the specific geometry. Rather than considering these conditions 
separately, we shall define the subclass of operators which satisfy evenness conditions at both 
faces. 

Definition 3.5. Let Ka be an element of the heat calculus ^f'^jf (M). Then Ka is an element 
of the even subcalculus if the following conditions hold: 
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i) At the face td, and in terms of the projective coordinates (T, H, £) (which are valid in 
the interior of this face), suppose that 

CO 

K A ^T- m J2^ d (^OT'; 

3=0 

then we require that 

4 d (-~,a = (-l)^ d (H,0 

for all j > 0. 

ii) At the front face S, and in the projective coordinates (p,£,£,w,y,z,z) from (I3.ip . 
valid in the interior of that face, suppose that 

CO 

K A ~ p- b - 3 J2 I w, y, z, z)p k ; 

k=0 

then we require that 

i4(^^-w,y,z,z) = {-l) k Kf{£,£,w,y,z,z). 

We shall denote the set of all operators which satisfy these conditions by ^ e !^if evn (-^0- 

There are a few different motivations for this definition. First, a straightforward calculation, 
see [221 §7-1] for the computation at td, shows that if Ka G C°°(M X M X [0, oo)), then its lift 
to M% satisfies both conditions i) and ii). Second, it is also not hard to check that the initial 
parametrix in the construction above can be chosen to satisfy both of these conditions. 
On the other hand, it is not immediately obvious that ^/*_ h . evn is actually a subalgebra, i.e. 
closed under composition. The remainder of this section is devoted to a proof of that fact, 
through the following two propositions. 

Proposition 3.6. Let P be any differential operator on M + x M x M which is a locally finite 
combination of the vector fields tdt, x 2 dt, xd x , x d y and d z with coefficients smooth in x 2 . 
Let P = P*P be its lift to M%. Then 

Po^_ h . evn (M)c*:_ h;evn (M). 

This is straightforward to verify in local coordinates. Note that this implies that if g is an 
even asymptotically admissible metric on M, then the higher order terms of tC preserve the 
even sub calculus. 

The more substantial result is the following: 

Proposition 3.7. 

c— n;cvn e— n;evn c— n;evn 

The technical proof is in the appendix, following the proof of composition for the full heat 
calculus. 

We now use these two propositions and follow through the parametrix construction of the 
previous section to conclude the following result. 

Proposition 3.8. The fundamental solution H for the heat equation associated to the 
Friedrichs extension of the (rescaled) Hodge Laplacian on a simple edge space (M,g) with 
even asymptotically admissible metric is an element of ^1'\. evn {M, e f\*). 
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4. Heat trace asymptotics 



Let us now turn to the form of the heat trace expansion, with particular attention to the 
implications of Proposition 13.81 on that expansion. 

First recall the diagonal Dq C M + x M 2 , and more importantly, its lift Dh to as 
pictured here. 



If 




D h C M 2 



td 



Figure 4. The diagonal hypersurface Dh in Ml. 



If ld : Dh M? is the natural embedding and the heat kernel for the Friedrichs extension 
of the Hodge Laplacian (which we consider as living on M? already, for convenience) , then by 
Proposition 13.41 t* k Hk is poly homogeneous on Dh, and hence so is tr t* k Hk, its pointwise trace. 
Our next goal is to obtain information about the terms which can appear in the expansion of 
this function. Let D t d = Dh D td, Dg = Dh D ff and D c = Dh D rf n If be the boundary faces 
of D h . 



Proposition 4.1. Let g be asymptotically admissible, 
index family is given by 



Then tvi* D H k G A% g {D h ) 



where the 



Gtd = -m + 2N , G s = -1 - b + N , G c = E u + E ri 



Moreover, if g is even, then Gs 



-1-6 + 2N, 



o- 



This Proposition is an immediate consequence of the polyhomogeneous structure of 
itself. The fact that the restriction of a polyhomogeneous distribution to the submanifold Dh 
is polyhomogeneous is straightforward and can be checked in local coordinates, but is also a 
special case of the 'pullback theorem' in [23J. Hence the form of the expansion in the general 
case follows directly from Proposition 13.41 except for the statement that the expansion at td 
has no odd terms. 

The vanishing of the odd terms in the expansion at td in general and at ff when g is even 
follows from Proposition 13.81 Indeed, the coefficients in the expansions of tiL* D Hk at each 
face of Dh are the pointwise traces of the coefficients in the asymptotic expansions of at 
the corresponding faces of M|. However, functions on td which are odd with respect to the 
reflection S 4 -S, and similarly functions on ff which are odd with respect to w h-> —w, 
must vanish on Dh, and thus have vanishing trace. It is a general fact, explained in [22j . that 
in the interior, always satisfies the evenness condition at td, while this evenness at ff only 
holds if g is even. 
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Our main interest, however, is not in the pointwise heat trace but in its integral, Tr H k . Note 
that .D/i can also be obtained as the parabolic blowup of M + x M at {t = x = x = 0, y = y} 
without any reference to the double space. Let 

vr c : D h M+ 

be the composition of the blowdown map Dh — > M + x M and the projection M + x M — > M + . 
Then Tr H k is the pushforward of ti i* D Hk by ir c : 

(4.1) TrH k (t) = (ir c ),(tTL* D H k )= f tr i* D H k (t,x,y, z) x~ f dV g . 

JMx{t} 

(The extra factor of x~f in the integrand is due to the rescaling $£..) The formula with the 
integral on the right is valid when t > since Dh is the same as the diagonal Dq of M + x M 2 
there, but that this expression must be interpreted carefully near t = because of the front 
face of D^. 

There is a good formalism for understanding pushforwards of polyhomogeneous distribu- 
tions, known as Melrose's Pushforward Theorem [23]. It states, roughly, that if / : X — )■ Y 
is a special type of map, known as a 6-fibration, between two manifolds with corners, and if 
u G Ap hg (X), then so long as the integrals along the fibres (or generalized fibres) make sense, 
we have that /*(«) G A^ hg (Y). Moreover, the index set of the image /*(«) is determined by 
the index set of u, and similarly the coefficients in the expansions of the image can be deter- 
mined in terms of (in general, regularizations of similar sorts of pushforwards) coefficients of 
u. For the reader's convenience, we review this theorem here in this specific geometric setting. 

We first note that the map tt c is indeed a 6-fibration. We refer to [23] for the precise 
definition of this condition (see also the appendix below), but do not repeat it here since it 
does not play a significant role. Using the two sets of projective coordinates on D^, arising 
from the coordinates (|3.ip and (|3,4p on M^, we find 

(4-2) <(t) = plpL 

It turns out that the formula for index sets of pushforwards is simpler if we write everything 
in terms of 6-densities. Thus suppose that is a density on Dh which is smooth up to all 
boundary faces and everywhere nonvanishing. A smooth 6-density p,b is, by definition, any 
density of the form (ptdPsPc)" 1 Po- There is a pleasant naturality of 6-densities with respect 
to blowups: a straightforward computation shows that 

Pb = Aj3^,((xt)~ 1 dtdxdydz) 

for some smooth nonvanishing function A. Note too that we can use pi, = /3p(t~ 1 x~f~ 1 dtdV g ) 
for this 6-density. 

Consider a polyhomogeneous 6-density K on L>^, so K = K/if, where K G A^ hg (Dh)- Now 
suppose that K has index family Q = (Gtd, Gg, G c ). In order that the integrals over the fibres 
of 7r c converge, one needs that all terms in the expansion of K at D c have strictly positive 
exponent, i.e. Hez > for all (z,p) G G c . Provided this integrability condition holds, then 
the pushforward theorem asserts that 

(w c ),K = (KMKp b ) = T(t)j, T G < hg (M + ), 

where the pushforward index set F = tt^(G) is defined as follows. If G t d = {( z j,Pj)} an d 
Gg = {wi,qi)}, then G( d = {(Zj/2,pj)}, G' s = {(wi/2,qe)} and F = G' td UG' s is the 'extended 
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union' of G' td and G' s 

(4.3) G' td UG' s = G' td UG' s U{((z,p + q + l):3(z,p)eG> td , and (z, q) G G' s }. 

We refer to [231 Eqn. (42)] for a proof. In our particular example we set K = 2p c ps (tr b^H^), 
so that T(t) = Tr H k (t); and obtain the heat trace asymptotics as above. 
We summarize all of this in the following central result. 

Theorem 4.2. Let (M m ,g) be a compact simple edge space, where g is an asymptotically 
admissible metric. Let H k denote the heat kernel for the Friedrichs extension for the Hodge 
Laplacian on k-forms. Then, as t — > 0, 

oo oo 

Tr H k (t) ^J2 A e t ^+22 Cit ~ + 22 Git ~ logt - 

1=0 1=0 l€3 

The index set 3 here is 

3 = {£ gN I £ + m-b even}. 

If the metric g is even, then 

00 oo 

Tr H k (t) ~ £ V* + f>tH + 

e=o t=o ley 

where 3' = if (m — b) is odd and 3' = No if {m — b) is even. 

We can say something about the coefficients which appear here. The coefficients bit of 
Pt d m+2i in the expansion of tr i* D H k (i.e. before integrating) are the standard local geometric 
quantities in the interior of M; on the other hand, the coefficients ai of Pg b ~ 1+i typically 
involve not only similar local quantities on the edge B and on all of M, but also global 
quantities over F. The coefficients Ap, Ci and Gi which appear above are integrals of these 
various coefficients which involve both local and global quantities. We do not comment on 
this further since it plays no role here, but intend to return to a closer examination of these 
coefficients in another paper. We refer to [29] for an example where global spectral invariants 
on F appear. 

We conclude this section by stating the implication of this main theorem on the meromor- 
phic structure of the zeta functions. Let P k denote the projector onto the nullspace of the 
Friedrichs extension of the Hodge Laplacian A/%, and define 

1 f°° 

Cfc(s) := tttt / t s - l Tr (H k - P k ) dt, Re (s) » 0. 

r (s) Jo 

By the usual arguments, Theorem 14.21 implies the following 

Proposition 4.3. Each Ck{ s ) extends meromorphically to the entire complex plane. Near 
s = 0, 

where 3 = {j G No | j + m — b even} . If g is also even, then 

OO j oo „ „ 

where 3' = if (m — b) is odd and 3' = No if (m — b) is even. 
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In particular, if g is even and m is odd, then Ck{ s ) * s regular at s = for every k. If g is 
even and both m and b are odd, then we even have Cfc(O) = 0. 



5. Analytic torsion and the metric anomaly 

We come, at least, to an examination of the analytic torsion of a compact simple edge 
space (M, g) with an asymptotically admissible edge metric. Recall that in terms of the zeta 
functions defined at the end of the last section, we define the analytic torsion zeta function 



Cat(s) :=-J2(-l) k kC(s,A h 

and then set 



2 

fc=0 



log T(M,g) = ^ 



Cat(s)- 

s=0 



In certain cases, the analytic torsion zeta function Cat may be regular at s = even though 
some or all of the individual zeta functions are not. For example, if m is even, then there 
is nothing in Proposition 14.31 which precludes £(s, A/%) from having a pole at s = 0. In the 
special case of isolated conical singularities (b = 0), it was shown in [13] that despite this, 
there is a cancellation of residues and Cat(s) is regular at s = 0. In fact, Dar's arguments 
also apply to more general situations where the conic cross-section F may have boundary or 
even be singular. A consequence of our work here, however, is that if F is smooth and m is 
odd, then each Cfe(s) is regular at s = 0, hence the analytic torsion is well-defined, without 
need for Dar's cancellation argument. 

Recall from the introduction that we consider a slightly different quantity than just the 
analytic torsion. The determinant line bundle over M, 



m 

detn*(M) := (g) (/\ m kerA 



^ ( _i)fc+i 
k ' 



fc=0 



has a natural Hermitian inner product and norm || • HdetWfiW) induced from the restriction 
to ker Afc of the L 2 norm on L 2 O fc . The more invariant object is the Quillen metric 



l(Af,g) := T ( M >9)\\ ■ ||detW(M)- 



We now study the variation of the analytic torsion on M as we vary g. Consider a family 
g^ of asymptotically admissible edge metrics on M, fi G (hq — e, no + e). Denote by Afe(//) 
the Hodge Laplacian on &-forms associated to g^. In order to make these operators act on 
the same Hilbert space, we conjugate by the natural local isometry 



= ^^J^. . L 2 n k (M,g^ — > L 2 n k (M, 9ll0 ). 



Thus consider the operators 

Afc(At) := 

as an operator on the fixed Hilbert space L 2 Q, k (M, g^). 



A k (fi) := T M o A k (fi) o T M 1 
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Revert now to the common notation Hfc(i) = e k and set H^it) = e k . The semigroup 
property of the heat kernel gives the identity 



e -tA k (p) _ e -tA k (p ) rt q I e -(t-s)A k (n ) e -sA k (n) 



Taking the limit fj, — > fj,Q yields 
d 



o ds y A* — Mo 

i 



(t- s )A fc ( M ) f Afc(^o) - A fc (/i) \ e - s A fc ( M ) ds _ 







M - Mo 



9^ 



Tr e " iA ^) 



Tr ( e -(*-s)MMo) CA fc (/io) J e- sAfc ^°) 



iTr A fc (/i ) e 



-tA fc (Mo) 



where the upper dot denotes the derivative in fj,. Evaluating Afc(/i) explicitly in terms of T^, 
we find 



d_ 



Tr f e -*A fe (M) 



t-Tr ff^oAfcOioJoe-^^or- 1 



t • Tr ( A fc (Ato) o T" 1 ^ o e -*A*(A«o)) _ t . Tr ( A*Oo) ° e 



-tA k ([i ) 



Using that bounded operators can be commuted under the trace, the second term on the right 
here becomes 

Tr (AM o f-% o e -*A*(MoA = Tr f e -*/2A fc (M„) Afc(M) y-l r e -*/2A fc (Mo) 



= Tr (f- % o A fc (Mo) o e~ tA *^ = _ Tr (V" 1 ^ o A fe (/i ) o e -* A ^o)) 
Consequently the first and the second terms cancel and we are left with 



tTr A fc (^o) ° e 



-tAfc( W ) 



If Pfe(jtx) is the orthogonal projection onto the kernel of Ajt(/u), then repeating the arguments 
in [26l pp. 152-153], we have 



£(_i)*. fc.Tr (e-^^-P,^)) 



(5.1) 



fc=0 



5 



fc=0 



Here * M is the Hodge star for ^ and := on forms of degree k. Put 



so that, by definition, logT(M,# M ) = ^| s=0 /(/•*, s )- 
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Now differentiate /(//, s) under the integral. Inserting (|5.1h . recalling that the heat trace 
decays exponentially and assuming Res 3> so that there are no boundary terms in the 
integration by parts, we have 



a 

(5.2 



"1 -1 /» 

r^ ( - 1)fc+1 r4r / ^ lTr ( c " tAfc04) - Pfe( ^)) *■ 



fc=0 

At this point we introduce the separate hypotheses on the family of metrics g^. First 
suppose that all the metrics are strongly asymptotically admissible, so that g^—g^ = 0{x b+1 ). 
Then, a* (a) = O (x 1+b ) as well. Using projective coordinates, and letting til ■ M + x M 2 — > M 

be the projection onto the left copy of M, we see that the lift (ttx ° /3)*ct^ is smooth on M| 
with 

(5-3) V 
Using the pushforward theorem once again, we deduce that 

oo oo oo 

J=0 j=l j=l 

On the other hand, suppose that each g^ is even in x. Then a^(x) is also even, and hence 
lies in the even subcalculus. Consequently, 

(oo 

{0 ' v ( -° \ 

By the pushforward theorem, we now have 

oo oo 

(5.5) Tr (aJe-* A ^)) ~ £ V"? + £ C^~l + £ G^~f log*, 

j=o i=o jea 

where 2J = if (m — 6) is odd and 3 = No if (m — b) is even. 

Under either of these two sets of assumptions, we have shown that for each k, 

Res s=0 j™ t s - lr Tr ( Q J ( e -' A *(") - P k (^))) dt = -Tr (a JP fc (//)). 

In view of the additional factor of s in (j5.2j) . we find 

(5.6) -logT(Jl^) = -£(-l) fc Ti: (^0*)) = 5- log II • Ude^wW 

^ fc=0 ^ 
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This proves our central result. 

Theorem 5.1. Let (M" 1 ,^) be a one-parameter family of asymptotically admissible edge 
metrics on a compact space with simple edge. Assume that either 

i) each metric is strongly asymptotically admissible 
with g^ = g^ + and |/i M | 9fl[) = 0(x l+h ), or 

ii) each metric is even asymptotically admissible. 

Then the family of analytic torsion norms \\ ■ ||^ g ^ is independent of the parameter \i: 

— II • ll^ 5 -n 

We anticipate that this metric independence will have some interesting applications. 

Appendix: The composition formula 

In this appendix we prove the main composition result for the incomplete edge heat calculus, 
and then show that composition preserves the even subcalculus. This is in some sense the 
most technically difficult part of the parametrix construction in the heat calculus. However, 
our method for understanding these compositions follows a now standard pattern in geometric 
microlocal analysis using Melrose's pushforward theorem for poly homogeneous distributions 
with respect to 6-fibrations between manifolds with corners. We begin by reviewing this 
general result and then describe how it implies the composition formula. The composition 
formula in the heat calculus is almost identical to the one in the conic setting proved by 
Mooers [24J, and is also very closely related to one for the pseudodifferential edge calculus in 
|21| ; we refer also to the survey [3] for more on these matters. Because of all these resources, 
we shall be brief here. 

To set the stage, however, suppose that A £ %>X{M) and B £ ^ c S h ' (M), with Schwartz 

kernels Ka and Kb, which are pushforwards under the blowdown (3 : M| — > M + x M 2 of 
poly homogeneous distributions ka and kb, respectively. Let C := A o B, with Schwartz 
kernel Kq lifting to kq on M?. Our main task here is to show that kq is poly homogeneous, 
with index sets at each face given as in the statement of the theorem in terms of those for ka 
and kb, respectively. 

To explain how this relates to the pushforward theorem, first observe that 

(5.7) K c (t,x,y,z,x',y',z')= f lK A {t-t',w,w')K B {t',w',w")dt'dV g {w') 

Jo Jm 

This can be rephrased as follows. Consider the triple-space x Ri", x M w x M w > x M w //, 
and the three projections 

tt c : M 3 x R+ x R+ M 2 x R t '+t", t", w, w',w") (w, w",t' + t"), 

(5.8) tt l : M 3 x R+ x R+ — » M 2 x R t ,, (t', t" , w, w' , w") (w, w',t'), 
7T R : M 3 x M+ x R+ -> M 2 x R t „, (t' , t" , w, w' , w") -> (w' , w" , t"). 

Assuming that all kernels are extended to vanish for negative times, and reinterpreting them 
as densities in a suitable way specified below, we can rewrite (|5.7p as 

K c = (ttc)* [kIKaitrKb) . 
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The main idea is that we define a heat triple-space M^; this will be obtained from (IR + ) 2 x 
(M) 3 by a sequence of blowups, and has the property that there are maps 

U L ,U C ,U R -Ml^Ml 

which 'cover' the three projections defined above. Lifting the composition to these spaces 
leads to the key formula 

(5.9) k c = (n c )* (n* L K A u* R K B ) . 

Thus it suffices to show that if ka and Kb are poly homogeneous, then their lifts to Mi are 
also poly homogeneous, and so too the product of these lifts, and most importantly, that the 
pushforward by He of this product is again polyhomogeneous on Mr. 



The pushforward theorem. We have already encountered a special case of this in the 
discussion preceding Theorem 14.21 First introduce some terminology. Let X and X' be two 
compact manifolds with corners, and let / : X — > X' be a smooth map. Let {Hi} and {Hj} 
be enumerations of the codimension one boundary faces of X and X', respectively, and let pi, 
p'j be global defining functions for Hi, resp. H'-. We say that the map / is a 6-map if f*pj is a 
smooth nonvanishing multiple of some product of nonnegative integer powers of the defining 
functions pi, or symbolically, 

fp'i = Aj n pf' j \ Aij > 0, e(i,j) G N U {0}. 

i 

This simply means that / respects the boundary structure of these two spaces, and in par- 
ticular maps each Hi into some Hj with constant normal order of vanishing along the entire 
face. 

The map / is called a 6-submersion if /* induces a surjective map between the 6-tangent 
bundles of X and X'. (The 6-tangent space at a point p of dX on a codimension k corner 
is spanned locally by the sections x\d xi , . . . ,Xkd Xk ,d y . , where the defining 

functions for the faces meeting at p and the yj are local coordinates on the corner through p.) 
If, in addition, the matrix e(i,j) defined above has the property that for each j there is at 
most one i such that e(i,j) ^ (this condition simply means that each hypersurface face Hi 
in X gets mapped into at most one H'- in X' , or in other words, no hypersurface in X gets 
mapped to a corner in X'), then / is called a 6-fibration. 

We already introduced the notion of a 6-density in §4. Let us fix smooth nonvanishing 
6-densities v\, on X and v' h on X' . 

Proposition 5.2 (The Pushforward Theorem). Let u be a polyhomogeneous function on X 
with index sets Ej the faces Hj of X. Suppose that each (z,p) £ Ej has Rez > if the index 
j satisfies e(i,j) = for all j (which means that Hj is mapped to the interior of X' ). Then 
the pushforward f*(uvb) is well-defined and equals hv' h where h is polyhomogeneous on X' and 
has an index family fb{£) given by an explicit formula in terms of the index family £ for X. 

Rather than giving the formula for the image index set in generality, let us describe it 
slightly informally but specifically enough for the present situation. If Hj x and Hj 2 are both 
mapped to a face H[, and if Hj 1 n Hj 2 = 0, then they contribute the index set Ej 1 + Ej 2 to 
H[. If they do intersect, however, then the contribution is the extended union Hj 1 UHj 2 . 



Rafe Mazzeo and Boris Vertman 



31 



The triple space. We now construct the reduced heat triple space M 3 h . The steps are dic- 
tated strictly by the requirements that the maps ttl, ttc, ttr all lift to 6-fibrations rl^, lie, I1r, 
and that the space M 3 h be as 'small' as possible. We revert to adapted boundary coordinates 
w=(x,y,z). _ 

The triple space M 3 h is a parabolic blowup of M 3 x R^ x RT,,, where the standard local 
coordinates in the singular edge neighborhoods of each of the three copies of M are (x, y, z), 
(x',y',z') and {x" , y", z"). First blow up 

F = {t' = t" = 0, x = x' = x" = 0, y = y' = y"}, 

parabolically with respect to both t' and t"; then blow up the resulting space [M 3 x R + x 
R + , F, {dt' ,dt"}] at the lift of O = {t' = t" = 0}; finally blow up the resulting space, parabol- 
ically in the respective time directions, at the lifts of each of the three submanifolds 

F C = {? = t" = 0,x = x" = 0,y = y"}, 

(5.10) F L = {t" = 0,x' = x" = 0,y' = y"}, 

F R = {t' = 0,x = x' = 0,y = y'}. 

Thus altogether, 

M 3 x R + x M+, F, {dt'}, {dt"}; O; F c , {dt', dt"}; F L , {dt"};F R , {dt'} . 

This is a rather difficult space to visualize, but one may 'see' part of it by ignoring the time 
directions; the spatial part of M 3 h is then exactly the same as the triple space appearing in 
the elliptic theory of edge operators, see [21 j . which can be pictured as in Figure 5. 





Figure 5. The spatial component of the triple space M 3 h . 

Here, (101), (Oil) and (110) label the boundary faces created by blowing up Fc,Fl and 
Fr, respectively. The face (111) is the front face introduced by blowing up F. We denote the 
defining function for the face (ijk) by Pijk- 

Now recall the projections ire, ttl arid itr defined in (|5.8|) . These induce projections lie, 
U L and Il R from M 3 h to the reduced heat space M r 2 h := [M 2 x R+; FD; {dt}] introduced in §3 
(see Figure 2). It is not hard to check that the choice of submanifolds that have been blown 
up ensures that these are in fact 6-fibrations. 

Denote the defining functions for the right, front and left faces of each copy of M r 2 h by 
{pio, Pu, Poi}> respectively. These lift via the projections according to the following rules 

^h(Pij) = PiOjPilj, 
(5.11) ^l(Pij) = PijOPijl, 

^R(Pij) = POijPlij- 
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Now consider the behaviour in the time variables. Let (3^ : M 3 h — > (M + ) 2 x M 3 be the 
blowdown map, tq the defining function for the face obtained by blowing up 0, and t', t" the 
defining functions for the two boundary faces in M r 3 h corresponding to {t' = 0} and {t" = 0}. 
Then we find 

((3^)*(t') 2 - 2 - 2 



(5.12) (/3 (3) )*(t") =t"~ 2 2 2 

(/3(3))*(t' + t") - - 2 



T r OPlllPnoPl01' 

: t"toPiiiPoiiPioi> 

T OPinPioi- 

Let /3 (2) : M 2 h -> M+ x M 2 be the blowdown map for the reduced heat space. The lifts 
(/3( 2 ))* (f), (/3( 2 ))* (i") and (/3( 2 ))* (i' + i") to n L (M 3 J, E^M* ) and n c (M 3 h ), respectively, 
are equal to Tp| 1; where T is the defining function for tf in M 2 h . Note 

n^°(/3 (2) )* = (/3 (3) )*o^,L,«- 

Consequently, in view of (|5.12p and (|5.1ip . we conclude 

n^(T) = r , 

(5.13) Ul(T) = r'rop 2 wl , 

U* R (T)=r"ropiof 

Using these data, we now derive the anticipated composition formula. Consider two 
elements of the heat calculus A £ *J?jf lf ' £rf (M) and 5 G ^'^'^'^ (M) . We con- 
sider the Schwartz kernel of each as a 'right density', -Ka(*' ,x,y, z,x' ,y' , z')dt'dx'dy'dz' and 
K B (t",x',y',z',x",y",z")dt"dx"dy"dz". Then their product on (M+) 2 x M 3 is 

KA(t', x, y, z, x , y', z')KB(t", x', y', z' , x", y" , z")dt' dt" dx' dy' dz dx" dy" dz" . 

The integral over dx 1 dy' dz' and over t' + t" = t gives Kaob(1-> x,y, z, x" ,y" , z")dt dx" dy" dz" . 
To put this into the same form required in the pushforward theorem, multiply this expression 
by dxdydz. 

Blowing up a submanifold of codimension n amounts in local coordinates to introducing 
polar coordinates, so that the coordinate transformation of a density leads to (n — l) st power 
of the radial function, which is the defining function of the corresponding front face. Hence 
we compute the lift 

(f3 {3) )*(dt'dt"dxdydzdx'dy'dz'dx"dy"dz") 

_ „2b+6.6+5 .6+3 .&+3_ ,.(3) _ n 2b+6 b+5 b+3 b+3 2 _/ // m„ \ „( 3 ) 

-Pin Pioi Pno Poii T o v y > - Pm PioiPiioPon r o r r {UPiik)v b , 
where is a density on Af 3 h , smooth up to all boundary faces and everywhere nonvanishing; 

i/j is a 6-density, obtained from by dividing by a product of all defining functions on M 3 h ; 
and (Hpijk) is a product over all (ijk) E {0, l} 3 . Furthermore, the infinite order vanishing 
of ka and Kb at T = implies that the product of the lifts vanishes to infinite order in 
Pioitot't" . Altogether, we obtain that 

(U* l k a ) (n* R K B ) (p (3) )*(dt'dt"dxdydzdx'dy'dz'dx"dy"dz") = p e +(' +1 (Up ljk ) Gvf \ 

where G is a polyhomogeneous function on M 3 h , vanishing to infinite order in Pioitot't" , with 
index sets E\f, E T f and E\{ + E T { at the faces (001), (100) and (010), respectively. Moreover, 
G has index sets E\f + £' and E' r{ + 1 at the faces (011) and (110), respectively. 
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Note that since ka does not vanish to infinite order on td, the lift I\* l ka is not polyhomoge- 
neous on MA. Fortunately, the other factor kb does vanish to infinite order there, and hence 
the product H* l ka • ITrKb is indeed polyhomogeneous on M r 3 h . Applying the Pushforward 
Theorem now gives 

(n c ), ((n^ A ) (n* R K B ) (^ 3) y(dt'dt"dxdydzdx'dy'dz'dx"dy"dz"^j 
(5-15) = (p^Y (K AoB (t,x,y,z,x",y",z")dtdxdydzdx"dy"dz") 

where is a 6-density on M 2 h , and (3^ : M 2 h — > M + x M 2 the corresponding blowdown 
map, and Q is a polyhomogeneous function on M 2 h , which vanishes to infinite order in T, has 
leading order t + £' in its asymptotic behaviour at ff , and with the index sets (P\{ , P r f) at the 
left and right boundary faces arising as extended unions (recall (|4.3p ) 

, . P u = E[,U(E u + £'), 

5-16 -, , 

Prf = £rfU(£ rf + £), 

By an argument similar to (|5.14p . we compute 

(5.17) (/9 (2) )* {dtdxdydzdx"dy"dz") = p\^ (pioPnPoiT) . 

Consequently, combining (|5.15p and (|5.17p . we deduce that {^)*K A oB = k-AoB has index 
sets P\{ and P r f at the faces (01) and (10) of M r 2 h , is p- b ~3+ f -+ f - times a smooth function at 
ff, and vanishes to infinite order in T. This proves 

Theorem 5.3. For index sets En and E' r{ such that En + E' r{ > —1, we have 
where 

P l t = E[ i U{E ]i +e), 
Prf = E li U{E' ri + £). 

The even subcalculus 

We now give the proof of Proposition lBTTl and show that composition preserves the subspace 
of operators which satisfy the parity conditions from Definition ^. 51 In fact, we refer to [22] for 
a proof that operators which satisfy these conditions near td are closed under composition, 
since the proof is exactly the same here. Thus we focus on behaviour near ff. It clearly 
suffices to consider operators with Schwartz kernels supported near some compact subset of 
the interior of the front face. Hence we are free to write out and examine the composition 
using projective coordinates. 

Note first that the pushforward by lie does not introduce logarithmic terms in the front face 
expansion of kaoB, since the kernel on M^ h is vanishing to infinite order at (101). Hence, for a a 
and kb with integer exponents in their front face expansions, k Ao b has a front face expansion 
with integer exponents as well, and we may check its evenness according to Definition 13.51 

As in the 'full' composition formula, we lift the integrand in the formula 

KAoB(t,w) = / / KA(t — t' ,w,w')KB(t' ,w' ,w) dt'dw' 
Jo Jm 
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to M^ h . For the left factor, use the coordinates (13, lj) . while on the right we introduce a further 
singular coordinate change 

y'-y" t> x> 

v = , cr = t>? = — > 

P t p 

where p = Ji. Following the substitutions in formula (18) of |16j . we obtain 

KAoB(p,t;,€,w,y,z,z) = 

1+5 f 1 f , p. 6 £' io — w „ . 

• ! 1 'kaUVI-o-, , ,-=^=,y + vp) 

\/l — a \/l — a \/l — a 



kb(px/&, "4=j "7=) ~~F="iy) d£ dv dz' . 
Jo Jo Jo 



Now expand the two factors to get: 

«a(PV1 -o", 7= = , 7= , 7= ,y + vp) 

s/l — o \/l — o \/l — o 

~ (pVT^)- 6 - 3 ^ ^(-=1=, -=!=, -^L, y + i,p)(pVT=7) fc , 
" Vl - cr V 1 - cr V 1 - cr 

kbo^, 4=. 4=' -?=' ~ o^r 6 - 3 E 5fc (4=> -t=> -7=' 

Jo Jo Jo ^ Jo Jo Jo 

By assumption, and l?2fc are even while ^fc+i an d B2k+i are odd in their third slots. 
Since each is smooth in the final variable, we expand further to obtain 

^fc( t= , 7= ,-==,y + vp) ~ > A( jfej -)(-==, -7==, 7= ,y,v)pr>. 

\/\-0 Jl-O JX-O ^ Vl - cr V 1 - V 1 - cr 

Note that clearly, 

A f j ? W ~ V ~ \ l 1\3 A I j ^ W-V ~ . 

7r^' 7T=P 7T=P y ' " v) = ^Tr^P 7r^' 7r^' y ' u) " 

Now define 

A [i\ := E A (kj)- 

k+j=i 

Then the two parity conditions together give that 

Inserting these expansions into the integral leads to an expansion of kaoB- Since the 
rescaled volume form x~^dV g is a smooth function of pv and p 2 , it does not affect these parity 
considerations. 

Collecting all the factors involving only o into a single function f(o), we get 



KAoB(p,S,S,w,y,z,z) -^ 3 V/ +t / / f((T)A [il (-^=, J. ,^==L,y,v) 

^ Jo J Jl-o Jl-o Jl-o 

xB k (-¥=,-^,^L,y)dodZ'dvdz'=-. p- b -' s V p* +k (A o B) itk (£, I w, y, z, \ 
Jo Jo Jo 
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Finally, observe that 

(Ao B) i)k (€,€,-w,y,z,z) 

1 f tr \a i £ £' w + v ~ \n i £ £ 



o 



\f \ - a' \J\- a' Vl - c' ' v 7 ^' V^' v 7 ^ 

= (-l)«-*(AoS) (i »({,{,«,,f,z,5). 
Hence the Schwartz kernel of A o £> satisfies the parity condition too. 
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